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We study anomalous scaling and multiscaling of two-dimensional crack profiles in the random fuse 
model using both periodic and open boundary conditions. Our large scale and extensively sampled 
numerical results reveal the importance of crack branching and coalescence of microcracks, which 
induce jumps in the solid-on-solid crack profiles. Removal of overhangs (jumps) in the crack profiles 
eliminates the multiscaling observed in earlier studies and reduces anomalous scaling. We find that 
the probability density distribution p{Ah{£)) of the height differences Ah{l) = [h{x -f i?) — h{x)] of 
the crack profile obtained after removing the jumps in the profiles has the scaling form p{Ah{(.)) = 

( {Ah^\e)y/^ ) ' ^^'^ follows a Gaussian distribution even for small bin sizes £. The 

anomalous scaling can be summarized with the scaling relation L2,^2(2/2))'-/^ I 
where {Ah^{L/2)y^^ ~ L«. 
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I. INTRODUCTION 

For over two decades, scaling of fracture surfaces has 
been a well studied, yet a controversial issue [H, |3| • Ex- 
periments on several materials under different loading 
conditions have shown that the fracture surface is self- 
affine Q , which implies that if the in plane length scales 
of a fracture surface are scaled by a factor A then the 
out of plane length scales (height) of the fracture sur- 
face scales by A*", where C is the roughness exponent. 
Many experiments on several materials including metals 
glass @, rocks 0] and ceramics Q have tested the 
scaling in three dimensions. The scaling regime in some 
cases has been quite impressive, spanning five decades 
in metallic alloys [8]. It is an interesting question as 
to whether the self-affinity measured so often can be re- 
placed by more complicated scenarios and whether in any 
particular setup and geometry the exponents are univer- 
sal as in the line-depinning scenario [1, ■ 

In two dimensions, recent studies have debated the 
picture of simple self-afhnity. In other words, for a 
two-dimensional crack profile h{x) one can look at vari- 
ous statistical measures including the dependence of the 
roughness on sample/system size and the scaling of var- 
ious moments of h(x). There is a discussion on whether 
the two-dimensional fracture surfaces would exhibit self- 
aflnne or multi-afhne scaling [H [H, [H, [11. Ref. [H 
argues that a crack line h(x) in two-dimensions is not 
self-affine; instead, it exhibits a much complicated multi- 
affine (or multiscaling) structure, with a non-constant 
scaling exponent for the g-th order correlation function 



C9W = (|/i(a;+^)-/i(x)|«)i/« - i^-i. In analogy to kinetic 
roughening of surfaces, it has been arg ued that fracture 
surfaces exhibit anomalous scaling [l5|: the global expo- 
nent describing the scaling of the crack width with the 
sample size is larger than the local exponent measured on 
a single sample [16, 17] . This means that the typical slope 
of h{x) develops an algebraic dependence on the system 
size L, and it is necessary to introduce two roughness 
exponents a global one {() and a local one (Ooc) whose 
difference measures the L-dependent extra lengthscale. 

In two dimensions, the available experimental results, 
mainly obtained for paper sampleSj_indicatc a roughness 
exponent in the range C^0.6-0.7P,i3,ll8|,ll9ll20|,il|. 
However, one should note that apparently one can mea- 
sure for various ordinary, industrial papers values that 
are significantly higher than ( = 0.7 The reasons 

for these discrepancies are not clear. It has also been 
noted that the roughness exponent is dependent on the 
crack velocity: at the onset of fast crack propagation the 
exponent makes a small jump from its value when the 
crack still grows in a stable fashion [23| . 

The theoretical understanding of the origin and uni- 
versality of crack surface roughness is often investigated 
by discrete lattice (fuse, central-force, and beam) mod- 
els. In these models the elastic medium is described by 
a network of discrete elements such as fuses, springs and 
beams with random failure thresholds. In the simplest 
approximation of a scalar displacement, one recovers the 
random fuse model (RFM) where a lattice of fuses with 
random threshold are subject to an increasing external 
voltage [23|. Using two-dimensional RFM, the estimated 
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crack surface roughness exponents are: C = 0.7 ± 0.07 
dl, Cioc = 2/3 [H, and C = 0.74 ± 0.02 [H]. Recently, 
using large system sizes (up to L = 1024) with exten- 
sive sample averaging, we found that the crack rough- 
ness exhibits anomalous scaling ^] . The local and global 
roughness exponents estimated using two different lattice 
topologies are: C/oc = 0.72 ± 0.02 and ( = 0.84 ± 0.03. 
Anomalous scaling has been noted in the 3D numerical 
simulations as well [2§]. The origins of anomalous scal- 
ing of fracture surfaces is not yet clear although recent 
studies [nl, [H, [l^ suggest that the origin of multiscaling 
and anomalous scaling in numerical simulations may be 
due to the existence of overhangs (jumps) in the crack 
profile. 

In this paper, we further quantify the influence of these 
overhangs in the crack profiles on multi-afflnc scaling and 
anomalous scaling of crack roughness exponents. In par- 
ticular, the questions we would like to address in this ar- 
ticle are the following: (i) whether anomalous scaling of 
roughness observed in numerical simulations is a result of 
these overhangs (or jumps) in the crack profiles, and (ii) 
whether removing the jumps in the crack profiles com- 
pletely eliminates multiscaling. This should then imply 
a constant scaling exponent Qoc such that the q-th order 
correlation function Cq{£) = {\h{x + i) - h{x)\iy^'^ - 
£'''°'=. It should be noted that Gaussian distribution for 
piAh{i)) has been noted in Refs. [HI, El IH, Ho] only 
above a characteristic scale where self-afRne scaling of 
crack surfaces is observed. In this study, we would like 
to further investigate whether removing these jumps in 
the crack profiles extends the validity of Gaussian prob- 
ability density distribution p{Ah{£)) of the height differ- 
ences Ah{£) = [h{x + £) — h{x)] of the crack profile to 
even smaller window sizes £. We also discuss the cases of 
open (OBC) and periodic boundary conditions (PBC), 
since the presence of the former might have an effect 
on whether "anomalous scaling" exists. For the PBC 
case, we show that the crack profiles can be collapsed to 
a "semi-circle law", a scaling ansatz followed by many 
stochastic processes that return to the origin [30|. The 
rest of the article consists of three sections: first we intro- 
duce the numerical details. In Section III, we go through 
all the numerical results, and finally Section IV presents 
the conclusions. 



II. MODEL 

We consider numerical simulations using two- 
dimensional random fuse model (RFM), where a lattice 
of fuses with random threshold are subject to an increas- 
ing external voltage . The lattice system we consider 
is a triangular lattice of linear size L with a central notch 
of length ao (unnotched specimens imply oq — 0). All 
of the lattice bonds have the same conductance, but the 
bond breaking thresholds, t, are randomly distributed 
based on a thresholds probability distribution, p{t). 
The burning of a fuse occurs irreversibly, whenever 



the electrical current in the fuse exceeds the breaking 
threshold current value, t, of the fuse. Periodic boundary 
conditions are imposed in the horizontal directions (x 
direction) to simulate an infinite system and a constant 
voltage difference, V, is applied between the top and the 
bottom of the lattice system bus bars. 

A power-law thresholds distribution p{t) is used by as- 
signing t — , where X e [0, 1] is a uniform random 
variable with density px^X) — 1 and D represents a 
quantitative measure of disorder. The larger D is, the 
stronger the disorder. This results in t values between 
and 1, with a cumulative distribution P{t) = t^/^. The 
average breaking threshold is < t >— 1/{D + 1), and 
the probability that a fuse will have breaking thresh- 
old less than the average breaking threshold < t > is 
P(< t >) = {1/{D + l)y^°. That is, the larger the D 
is, the smaller the average breaking threshold and the 
larger the probability that a randomly chosen bond will 
have breaking threshold smaller than the average break- 
ing threshold. 

Numerically, a unit voltage difference, V = 1, is set be- 
tween the bus bars (in the y direction) and the Kirchhoff 
equations are solved to determine the current flowing in 
each of the fuses. Subsequently, for each fuse j, the ra- 
tio between the current ij and the breaking threshold tj 
is evaluated, and the bond jc having the largest value, 
maxj^, is irreversibly removed (burnt). The current is 
redistributed instantaneously after a fuse is burnt imply- 
ing that the current relaxation in the lattice system is 
much faster than the breaking of a fuse. Each time a 
fuse is burnt, it is necessary to re-calculate the current 
redistribution in the lattice to determine the subsequent 
breaking of a bond. The process of breaking of a bond, 
one at a time, is repeated until the lattice system falls 
apart. 

Using the algorithm proposed in Ref. [IH, we 
have performed numerical simulation of fracture up 
to system sizes L = 512 for unnotched samples and 
up to L = 320 for notched samples. Our simula- 
tions cover an extensive parametric space of {L, D 
and ao) given by: L = {64,128,192,256,320,512}; 
D = {0.3,0.4,0.5,0.6,0.75,1.0}; and gq/L = 
{0,1/32,1/16,3/32,1/8,3/16,1/4,5/16,3/8}. A mini- 
mum of 200 realizations have been performed for each 
case, but for many cases 2000 realizations have been used 
to reduce the statistical error. 



III. CRACK ROUGHNESS 

A. Crack width 

Once the sample has failed, we identify the final crack, 
which typically displays dangling ends and overhangs 
(see Fig. [iJ. We remove them and obtain a single val- 
ued crack line h^, where the values of x G [0, L]. For 
self-affine cracks, the local width, w{l) = {J2xi^x — 
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(1/0 hxYY^"^, where the sums are restricted to re- 
gions of length I and the average is over different real- 
izations, scales as w{l) ~ K for / ^ L and saturates to 
a value W — w{L) ^ L'^ corresponding to the global 
width. The power spectrum S{k) = {h^h-k) / L, where 
hk = J2x^xexpi{2Trxk/L), decays as S(k) ~ 
When anomalous scaling is present [isl . fig . [l7| , the ex- 
ponent describing the system size dependence of the sur- 
face differs from the local exponent measured for a fixed 
system size L. In particular, the local width scales as 
w{£) ~ iCioaj^C-Cioc ^ gQ ^j^g global roughness W scales 
as L'^ with ( > Qioc- Consequently, the power spectrum 
scales as S{k) - fc-(2Cioc-i-i)7,2(c-C,oc). 
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FIG. 1: (Color online) A typical crack in a fuse lattice system 
of size L X L with L = 512. This crack, identified as (a) in 
the figure has dangling ends, which are removed to obtain a 
single valued crack profile h{x), identified as (b) in the figure. 
This final crack h{x) possesses finite jumps that arise due to 
the solid-on-solid projection to obtain a single- valued fracture 
surface. The inset shows a zoomed portion of the crack. 

Figured^ presents the scaling of local and global crack 
widths in systems with different disorder values and an 
initial relative notch size of aq/L = 1/16. The slopes 
of the curves presented in Fig. [5^ suggest that a lo- 
cal roughness exponent Czoc = 0.71 that is independent 
of the disorder. The global roughness exponent is esti- 
mated to be C = 0.87, and differs considerably from the 
local roughness exponent C;oc- The collapse of the data 
in Fig. [2j) clearly demonstrates that crack widths fol- 
low such an anomalous scaling law. Notice that we have 
scaled away the amplitudes of the roughness for all the 
different D to achieve the maximal data collapse to illus- 
trate the universality. In the range of D considered here 
the amplitudes vary by about 20 %. The inset in Fig. 

reports the data collapse of the power spectra based 
on anomalous scaling for different disorder values. This 
collapse of the data once again suggests that local rough- 
ness is independent of disorder. A fit of the power law 
decay of the spectrum yields a local roughness exponent 
of Cioc = 0-74. This result is in close agreement with the 



real space estimate and we can attribute the differences 
to the bias associated to the methods employed (32| . 
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FIG. 2: (Color online) (a) Scaling of local and global widths 
w{l) and W of the crack for different system sizes L — 
{64, 128, 256, 320}, disorder values D and a fixed ao/L = 1/16 
value (top). The local crack width exponent (^loc ~ 0.71 
is independent of disorder and difi'ers considerably from the 
global crack width exponent C, — 0.87. (b) Collapse of the 
crack width data using the anomalous scaling law (bottom). 
Lc — {L — ao) is the effective length of the crack profile. Col- 
lapse of the data including a disorder dependent prefactor 
A{D) for a given disorder value implies that local and global 
roughness exponents are independent of disorder. The inset 
shows collapse of power spectrum S{k) using the anomalous 
scaling law with i^joc = 0.71 and C, = 0.87. The slope in the 
inset defines the local exponent via —{2(^ioc 4- 1) = —2.48. 
(a)-(b) present a total of 20 data sets. 



B. Anomalous Scaling 

The scaling properties of the crack profiles h{x) can 
also be studied using the probability density distribution 
p{Ah{i)) of the height differences A/i(£) = [h{x + i) - 
h{x)] of the crack profile between any two points on the 
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reference line (x-axis) separated by a distance (.. Assum- 
ing the self-affine property of the crack profiles implies 
that the probability density distribution p{Ah{£)) follows 
the relation 



piAh{e)) ^ (A/i2(£))-i/2 / 



Ah{e) 



(A/l2(^)}l/2 



(1) 



where {Ah'^{f)y^'^ denotes the width of the height differ- 
ence Ah{£) over a length scale i. 

Since for PBC the periodicity in crack profiles is anal- 
ogous to return-to-origin excursions arising in stochastic 
processes, we propose the following ansatz for the local 
width (A/i2(£))i/2 in height differences Ah{£) 



{Ah^i)) 



1/2 _ 



(A/i2(L/2))V2^(_ 



(2) 



The function (j) ( | is 



with (A/i2(i/2))i/2 = L(_ 

symmetric about £ = L/2 and is constrained such that 

(j) (^j^^ = at ^ = and £ = L, and (/) (^-^) = 1 at 

£ = L/2. Based on these conditions, a scaling ansatz of 
the form 



(A/i2(L/2))i/2 



1/0, 



ii-L/2y 
{L/2Y 



= 1 



(3) 



similar to stochastic excursions or bridges can be pro- 
posed for (A/i2(£)^i/2^ which implies a functional form 



L/2 



L/2)\ 



L/2 ) 



(4) 



for (f> (172) that is satisfied to a good approximation by 
our numerical results. This scaling ansatz implies anoma- 
lous scaling when 7^ C- Upon further simplication, 
Eq. (g]) results in 



(5) 



which along with (A/i2(L/2))i/2 ^ ^nd Eq. ^ il- 
lustrates how anomalous scaling appears in the scaling 
of local widths (Ah^{£)y/^ , and how local and global 
roughness exponents Qoc and C can be computed based 
on numerical results. 

Figure [3] presents the scaling of (Ah'^{£)y^'^ based 
on the above ansatz (Eq. ([3])). The collapse of 
the (A/i2(^)}i/2/(A/i2(L/2))i/2 data for different system 
sizes L and window sizes £ onto a scaling form given 
by Eq. ^ with C/oc — 0.64 can be clearly seen in Figs. 
[3](a)-(c). In particular. Figs. [2Ia)-(c) present the data for 
unnotched and notched samples with varying amounts of 
disorder D and relative crack sizes qq/L. The collapse of 
the data for varying amounts of disorder (0.3 < D < 1) 
and relative crack sizes (0 < ao/L < 3/8) can be clearly 



seen in these figures and suggests that local roughness ex- 
ponent Qoc is independent of disorder, at least for the dis- 
order ranges considered here. It is an interesting question 
as to why the Qoc from collapsing the average crack pro- 
files does not agree with the value from the local width, 
demonstrated in Fig. O 

Figures |li;a)-(b) present the scaling of (A/i2(L/2))i/2 
for various notched and unnotched samples with varying 
amounts of disorder and relative crack sizes. The data 
presented in these figures shows that {Ah'^{L/2)y/^ - 
with ^ — 0.87 in agreement with the previously given 
value for the global width exponent. In Figure IDJb) one 
can note that there is a ao/i-dependent amplitude and 
the data follow the 0.87-exponent at fixed uq/L. Since 
there exists a significant difference between the global 
and local roughness exponents (C and Qoc, respectively), 
here again we can conclude that crack profiles obtained 
using the fuse models exhibit anomalous roughness scal- 
ing. 

The questions that we would like to resolve in the fol- 
lowing are whether this anomalous scaling and multi- 
affine scaling of crack surface roughness are a conse- 
quence of the jumps in the crack profiles induced by 
the crack overhangs (see Fig. [T]). As shown in Fig. [51 
removal of jumps from an initially periodic crack pro- 
file h(x) makes the resulting crack profile hNp(x) non- 
periodic, where the subscript NP refers to nonperiod- 
icity of the profiles. A direct evaluation of the rough- 
ness exponent using these nonperiodic profiles can be 
made. However, such an evaluation of roughness exhibits 
finite size effects for window sizes £ > L/2. Alterna- 
tively, the roughness of these resulting nonperiodic pro- 
files can be evaluated by first subtracting a linear profile 



hlin{x) 



INP 



(0) 



{hNp{L)-hNp{0)) , 



from the nonpe- 



riodic profile ftjv p (x) , and then evaluating the roughness 
of the resulting periodic profile hp{x). In the following, 
we consider the scaling of hp{x). 

Figure [S] presents the scaling of crack width w{£) with 
window size £ for crack profiles without the jumps. The 
data presented in Fig. [6^ suggests that local and global 
roughness exponents {Qoc — 0.74 and C, = 0.80) are not 
the same even after removing the jumps in the crack pro- 
files, although the difference between these exponents is 
small. We have also investigated the power spectra S{k) 
of the crack profiles without the jumps in the crack pro- 
files (see Fig. [5}d). An excellent collapse of the data is 
obtained using the anomalous scaling law for power spec- 
trum with 2{C,—C,ioc) = 0.1 and Ooc = 0.74. This result is 
consistent with the exponents measured using the crack 
widths as in Fig. [S^. 

In addition to the above two methods (crack width 
scaling and power spectrum method) used for estimat- 
ing the local and global roughness exponents, we also 
used the scaling ansatz proposed in Eq. ([3]) to estimate 
the local and global roughness exponents. Since the dif- 
ference between the local and global exponents is small, 
alternate ways of measuring these exponents provide a 
sense of reliability into these estimates. Figures [71[a)-(c) 
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FIG. 3: (Color online) Scaling of {A/i^(^))^''^ with window size £. (a) Un-notched fuse simulations for L = {64, 128,256, 512} 
and D = 1.0 (top left), (b) Notched fuse simulations with various disorder values D — {0.3,0.4,0.5,0.6,0.75}, system sizes 
L — {64, 128, 192, 256, 320, 512}, and a fixed notch size ao/L — 1/16 (top middle). A total of eighteen data sets are plotted in 
figure, (c) Notched fuse simulations for various relative crack sizes ao/L = {1/32,2/32,3/32,4/32,6/32,8/32, 10/32, 12/32}, 
system sizes L = {64, 128, 256, 320}, and a fixed disorder of D = 0.6 (top right). A total of twenty four data sets are plotted 
in figure. 



present the scaling of crack profiles hp{x). The collapse 
of the (A/i|,(£))i/V(A/i|,(L/2))i/2 data for different sys- 
tem sizes L and window sizes £ onto a scaling form given 
by Eq. ([3]) with C/oc = 0.62 can be clearly seen in Fig. 
[7ja). This value is fairly close to the 0.64 quoted before. 
In addition, the collapse of the data presented in Fig. 
m.c) for {Ah'^p{i)y/y {Ah'^p{L/2)y/i demonstrates that 
multi-affine scaling of fracture surfaces arises because of 
overhangs (jumps) with certain statistics in the crack pro- 
file and removal of these jumps in the crack profiles com- 
pletely eliminates multiscaling of fracture surfaces. 

The simple scaling (no multiaffinity) is also evident 
through the scaling of (A/ip(L/2))^/^ presented in Fig. 
[7]Jb) . The slopes of the data for moments g = 1 to 6 of 
Ahp{L/2) are identical. An interesting observation to be 
made is that (A/i|,(L/2))i/« - with C = 0.80 whereas 
the local roughness exponent as obtained from Figs. [Tl^a) 
and (c) is Qoc = 0.62. A similar behavior is observed even 
when the linearity in the profile is not subtracted: the 
scaling of these nonperiodic profiles hpfp(x) is in agree- 
ment with that obtained for periodic profiles for window 
size £ < L/2 although finite size effects are observed when 
window sizes £ > L/2 are considered. The difference in 
these exponents even after removing the jumps caused by 
overhangs in the crack profile indicates that anomalous 
scaling is present in two-dimensional fracture simulations 
using the fuse models and this anomalous scaling is not 
due to the jumps in the crack profiles. 

C. The case of open boundaries 

It is interesting to compare the PBC case with that of 
open boundary conditions. Figure [5] presents the scal- 
ing of crack widths for fuse lattice simulations with open 
boundary conditions. The data in Fig. [5^ indicates that 
local roughness exponent is Qoc — 0.75. However, the 
data for different system sizes does not collapse, which 



is an indication of anomalous scaling. The inset in Fig. 
[5^ shows that a simple L-dependent shifting of the data 
achieves a perfect collapse of the data with possible finite 
size deviations for window sizes £ approaching the system 
size. Figure Hb presents the scaling of w{£) for crack pro- 
files without the jumps. Even after removing the jumps 
from the crack profiles, the crack widths data does not 
collapse onto a single curve. This suggests that removal 
of overhangs in the crack profile does not eliminate this 
apparent anomalous scaling of crack roughness even for 
open boundary conditions. 

On the other hand, removing the jumps in the crack 
profiles once again completely eliminates the multiscal- 
ing. Figure [5^ presents the scaling of g-th order corre- 
lation function Cq{£) = {\h{x + £) — h{x)\'^y^'' measured 
using the original crack profiles. Multiscaling below a 
characteristic length scale can be clearly seen in Fig. |9^. 
The data in Fig. [5Jd represents the scaling of q-th order 
correlation function Cq{£) measured after removing the 
jumps in the profiles. The Figure shows that the plots 
for different crack profile moments q are parallel to one 
another, and thus the removal of jumps in the crack pro- 
files eliminates multiscaling. A collapse of these plots is 
shown in the inset and the local roughness exponent is 
estimated to be C,ioc — 0.72, close to the PBC value. 

In the following, we finally investigate the probability 
density p{Ah{£)) of height differences Ah{£). In Refs. 
[Til . the p{Ah{£)) distribution is shown to follow a 
Gaussian distribution above a cutoff length scale and the 
deviations away from Gaussian distribution in the tails 
of the distribution have been attributed to finite jumps 
in the crack profiles. A self-affine scaling oi p{Ah{£)) as 
given by Eq. ([1]) implies that the cumulative distribution 
P{Ah{£)) scales as P{Ah{£)) ~ P(A/i(^)/(A/i2(£))i/2). 
Figure llOf a) presents the raw data of cumulative prob- 
ability distributions P{Ah{£)) of the height differences 
Ah{£) on a normal or Gaussian paper for bin sizes £ L. 
As observed in Refs. [ll|,[ll, and in Ref. 20] Fig. [TUIa) 
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FIG. 4: (Color online) Scaling of {Ah'^{L/2)y/^ with system 
size L. For notched samples, we use the efltective length of 
the crack profile Lc = L - oq. (a) Scaling of {Ah?{L/2))^^'^ 
is shown for unnotched samples and for samples with a fixed 
relative notch size of ao/L = 1/16 having varying amounts of 
disorder D. (b) Scaling of (Ah^ for samples with 
varying notch sizes and a fixed disorder of D = 0.6 (bottom). 



shows large deviations away from Gaussian distribution 
for these small bin sizes. However, for moderate bin sizes, 
the distribution is Gaussian with deviations in the tails 
of the distribution beyond the 3a = 3(A/i2(^))i/2 \[^[^ 
(data not shown in Figure). Removing the jumps in the 
crack profiles however collapses the P{Ahp{£)) distribu- 
tions onto a straight line (see Fig. [TUTb)) indicating the 
adequacy of Gaussian distribution even for small win- 
dow sizes £. Indeed, Fig. [TOTb') shows the collapse of the 
P{Ahp{£)) data for a system size L — 512 with a variety 
of bin sizes 2 < £ < L/2. Removing the jumps in the pro- 
files not only made the P{Ahp{£)) distributions Gaussian 
even for small window sizes £ but also extended the va- 
lidity of P{Ahp{£)) Gaussian distribution for moderate 
bin sizes to a 4cr = 4:{Ah%{£)y/^ (99.993% confidence) 
limit. 
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FIG. 5: (Color online) Figure shows a typical single valued 
crack profile h{x) with jumps based on solid-on-solid projec- 
tion scheme (identified as (a)). Removing the jumps in the 
crack profile h{x) makes it a non-periodic profile (identified 
as (b)). Subtracting a linear profile from this non-periodic 
profile results in a periodic profile (identified as (c)). 



IV. DISCUSSION 

In summary, we have here considered the nature of the 
roughness of the crack surfaces in the two-dimensional 
RFM. The results presented here indicate universality 
of local roughness exponent for both notched and un- 
notched samples with different disorders D in the range 
0.3 < D < 1.0 and for different relative crack sizes qq/L. 
This is true both for open and periodic boundary condi- 
tions. 

The results indicate that anomalous scaling of rough- 
ness is a generic feature of two-dimensional fracture in the 
fuse model. This is in contrast to e.g. the beam model 
[sst . where the global and local exponents are equal. The 
difference of the global and local exponents arises due to 
an additional lengthscale, which scales as a power-law 
of the system size L. We further investigated whether 
anomalous scaling of roughness is an artifact of presence 
of large jumps in the tails of p{Ah{£)) distribution. To 
do this, we considered the Ah{£) data that is only within 
±3(7 range of mean of p{Ah(£)) distribution, and com- 
puted the corresponding (A/i^(^))^/^ for various window 
sizes £. However, the data even from these truncated 
p{Ah{£)) distributions showed anomalous scaling. We re- 
peated our investigation with ±2cr range as well, but with 
a similar result. This indicates that anomalous scaling of 
roughness is not due to the tails of p{Ah{£)) distribution 
and persists in the mean as a function of L. 

Our results provide a concrete proof that the apparent 
multi-scaling of crack profiles observed in Ref. [lj| is an 
artifact of jumps in the crack profiles that are formed due 
to the solid-on-solid approximation used in extracting the 
crack profiles. The removal of these jumps from the crack 
profiles completely eliminates this apparent multi-scaling 




FIG. 6; (Color online) Scaling of crack profiles without the 
jumps, (a) Scaling of crack width w{£). The local and global 
roughness exponents appear to be different although the dif- 
ference is considerably smaller than that with the jumps in 
the profiles, (b) Scaling of power spectra of crack profiles. 
An excellent collapse of the data is obtained using anomalous 
scaling law with 2(^ — Cioc) =0.1 and C^ioc = 0.74. Removal of 
overhangs in the crack profiles does not appear to eliminate 
the anomalous scaling of crack roughness in fuse models. 



of crack profiles. Furthermore, removing these jumps in 
the crack profiles extends the vahdity of Gaussian prob- 
ability density distribution p{Ah{£)) of the height differ- 
ences to even smaller window sizes £ and to a range (of 
4cr = 4(A/ip(£))^/^) well beyond that observed in earlier 
studies. 

In conclusion, though the RFM is a "toy model" of 
(two-dimensional here) fracture, it still poses interest- 
ing issues and can be used to study questions that are 
also relevant for experiments. Our numerical results pre- 
sented here raise three basic theoretical questions related 
to the morphology of two-dimensional RFM fracture sur- 
faces that still remain to be answered. First, why does 
the extra lengthscale that leads to anomalous scaling 
have to be algebraic? Second, models explaining the 
dynamics in the final avalanche (unstable crack propa- 
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FIG. 7: (Color online) (a) Scaling of {Ah%{£)y^'^ with window 
size £ (top); (b) Scaling of (A/ip(L/2))^''^ with system size 
L (middle); (c) Scaling of (A/ip(^))^^'' with window size £ 
for different moments q of crack profiles without overhangs 
(bottom). 
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FIG. 8: (Color online) Scaling of crack width w{£) with win- 
dow size £ for open boundary conditions, (a) Scaling of w(^) 
for crack profiles with jumps. The inset presents the data 
shown in the main figure after a L-dependent shift is applied. 
A power law fit to the data estimates the local roughness ex- 
ponent to be = 0.75. (b) Scaling of 'w{£) for crack profiles 
obtained after removing the jumps, plots in figure (b) indicate 
that even with the removal of overhangs in the crack profile 
does not eliminate this apparent anomalous scaling of crack 
roughness. 



gation) and the roughness exponent values would be the- 
oretically interesting to develop. Third, how strong is 
the universality of the roughness exponents for disorders 
very different from the ones used here, in particular, to 
those leading to percolative damage or finite densities of 
infinitely strong fuses? Finally, in addition to present- 
ing results that explain the origins of apparent multiscal- 
ing and anomalous scaling, we have also made a connec- 
tion between periodic fracture surfaces and excursions of 
stochastic processes. 
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